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We derive a hydrodynamical description of the eigenvalues of the chiral Dirac spectrum in the 
vacuum and in the large N (volume) limit. The linearized hydrodynamics supports sound waves. 

The stochastic relaxation of the eigenvalues is captured by a hydrodynamical instanton conhguration 
which follows from a pertinent form of Euler equation. The relaxation from a phase of localized 
eigenvalues and unbroken chiral symmetry to a phase of de-localized eigenvalues and broken chiral 
symmetry occurs over a time set by the speed of sound. We show that the time is Ar = 7rp(0)/2/3A 
with p(0) the spectral density at zero virtuality and /? = 1, 2,4 for the three Dyson ensembles that 
characterize QCD with different quark representations in the ergodic regime. 

PACS numbers: 12.38Aw, 12.38Mh, Tl.lOPm 


1. Introduction. QCD with light quarks breaks 
spontaneously chiral symmetry. As a result, the light 
quarks transmute to massive constituents which make up 
for most of the visible mass in our universe. Empirical 
evidence for the spontaneous breaking of chiral stonmetry 
is in the form of light pions and kaons in nature [l| . Ded¬ 
icated first principle QCD lattice simulations with light 
quarks have established that chiral symmetry is spon¬ 
taneously broken with a finite chiral condensate and an 
octet of light mesons 0 . 

The spontaneous breaking of chiral symmetry is char¬ 
acterized by a large accumulation of eigenvalues of the 
Euclidean Dirac operator near zero virtuality Q. This 
phenomenon signals the onset of an ergodic regime in the 
chiral Dirac spectrum. Other regimes where the light 
quarks diffuse or undergo ballistic motion can also be 
identified Q. In many ways light quarks interacting via 
colored Yang-Mills fields behave like disordered electrons 
in metallic grains. 

The essentials of the ergodic regime are captured by 
a chiral random matrix model Q. In short, the QCD 
Dirac spectrum near zero virtuality only retains that the 
QCD Dirac operator is chiral or block off-diagonal, with 
random entries that are sampled from the three univer¬ 
sal Dyson ensembles Q with a Gaussian weight thanks 
to the central limit theorem. The spectrum corrections 
are also generic and follow from the neighboring diffusive 
regime with the light quark return probability falling like 
a power law in proper time Q. Both regimes are sepa¬ 
rated by the Thouless energy SQ- 
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In this letter we develop a hydrodynamical descrip¬ 
tion of the eigenvalues of the QCD Dirac operator in the 
diffusive regime much along the lines of our recent 
studies of the eigenvalues of the Polyakov line at large 
number of colors We will use this derivation to ob¬ 
tain the following new results: 1/ A hydrostatic solution 
for the spectral density beyond Wigner semi-circular dis¬ 
tribution; 2/ a hydrodynamical instanton that captures 
the stochastic relaxation of the eigenvalues of the QCD 
Dirac operator for low virtuality; 3/ a dynamical relax¬ 
ation time for restoring and/or breaking spontaneously 
chiral symmetry directly from the QCD Dirac spectrum; 
4/ an estimate of this time for the dyon liquid model. 

2. Chiral Dirac spectrum. The original chiral ran¬ 
dom matrix model partition function for the eigenvalues 
of the chiral Dirac spectrum was initially suggested as 
a null dynamical assumption for the generic analysis of 
the macroscopic chiral moments in the instanton liquid 
model However and more importantly, it was noted 
that this assumption provides a universal description of 
the chiral moments in the microscopic limit li,[i3- For 
QCD with Nf qiiarks of equal masses m in the complex 
representation [1,0 


r 

Z2[m] = / rfF n 

•> f=i 

Here T is a symmetric random complex matrix 

capturing the random hopping between N-left and N- 
right zero modes. o was generalized to all Dyson ensem¬ 
bles with /3 = 1,2,4 corresponding to quarks in different 
representations in [1|, 
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I- N N 

Zp[m] = / l[\X^, - X^/l[xnxl +mT^e-^^Ux,i2) 

i<j i 

with a = f3{v + 1) — 1. v measures the topological asym¬ 
metry of T, that is the difference between the number of 
zero modes and anti-zero modes for instance. Here Xf are 
the eigenvalues of the squared Wishart matrix W = T'^T. 
The overall normalization in are omitted. 

The Gaussian measure is generic at large N by the 
central limit theorem, with the parameter a fixed by the 
chiral condensate. In the microscopic limit whereby the 
Dirac spectrum near zero-virtuality is magnified so that 
NXi = 1, the interactions between the eigenvalues as me¬ 
diated by the gauge-helds are chaotic. As a result only 
the chiral structure of the Dirac Hamiltonian and the 
symmetry of the matrix entries under time-reversal are 
relevant [Hi- This is the universal regime of ergodic- 
ity shared by most disordered electronic systems in the 
mesoscopic limit [l^ . 

Alternatively, ([2]) can be regarded as the normalization 
of the squared and real many-body wave-function 

N 

Zp[m]= I l[dK\^o[X^]f (3) 

i—1 

which is the zero-mode solution to the Shrodinger equa¬ 
tion = 0 with the self-adjoint squared Hamiltonian 


Here pjj is the Hilbert transform of p 

[p\h = Ph{X) = -V j dX' (8) 

The canonical pair, 7r(A) and p{X) satisfies the Poisson 
bracket {7r(A), p{X')} = 6{X — X'). We identify the collec¬ 
tive fluid velocity with v = d\n and re-write (IHl) in the 
more familiar hydrodynamical form 

H Ki J (iAp(A) -f u[p]) Ri J d6p{X)\v + iAf (9) 

ignoring ultra-local terms. The equation of motion for p 
yields the current conservation law dtp = —2dx (pv), and 
the equation of motion for v gives the Euler equation 

dtv = {H, v} = —d\ -I- — d\A — Ad\lnp 

-|-7r/3[Ap]//(A)-|-7r/3[Ap]//(-A)] (10) 

All the relations hold for large but finite N, to allow for 
a smoothening of the density. 

4. Hydrostatic solution. The hydrostatic solution 
corresponds to the minimum of ([51) with r; = 0 and thus 
A(A) = 0. In terms of the chirally symmetric combina¬ 
tion p^{X) = p{X) + p{—X) we have 


N 

Hq = 'y ^ {—di -|- a^) {di + a.i) (4) 

i=l 

with di = d/dXi and the pure gauge potential a.^ = diS. 
The mass parameter is 1/2. The Vandermond contribu¬ 
tion A = ni<j induces a diverging 2-body part. It 
is convenient to re-absorb it using a similarity transfor¬ 
mation by defining = '^q/xTK and the new Hamitonian 

H = {1/VA)HoVA (5) 

3. Hydrodynamics. We can use the collective co¬ 
ordinate method in [l3l| to re-write in terms of the 
density of eigenvalues as a collective variable p(A) = 
~ ^i)- After some algebra we obtain 


H = J dXidxTTpdxir + pu[p]) 

with the potential-like contribution u[p] = A^ 


( 6 ) 


A= ^X-^[ph{X)-ph{-X)] 
+dxln^--^- 


2A A^ -I- 


( 7 ) 


0 = 


^A-^p/,(A) 


NfX 

X^+m^ 


( 11 ) 


In the large N limit only the terms in the first line survive 
with 


pI{X) = —J--X^ (12) 

TT V a 

which is Wigner Semi-circular distribution within the 
support |A| < 2lx/a. To correct it in 1/N we define 
pX ~ pX _|_ p_^ subject to the condition f pidX = 0, so 
that 


- d7rpij^(X) = j(a + 2Nf) - (13) 

The closed form solution is 


/\\ Q:-h2JVf |A| c-/\9 a! \ 

pi(A) =--— -S{X)-\ -p=J(A^-4/a) 


/Syi/o 


a + 2Nj — 


2^47^ 


jd'Kx/ 4/a — A^ 


(14) 
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with the support unchanged. We note that (|T4)) shows 
a delta-function accumulation at zero virtuality with the 
negative strength a + 27V/. 

The general solution to (HID for large but finite TV can 
be sought in the massless case or m = 0, by multiplying 
(ED by 2p^A(A^ and taking the Hilbert transform. 

The result is 

apN {\^p% - A^^ - ^A 

+\dxp% - (a + 2Nf)p% = 0 (15) 

To solve it consider the extension of the resolvent 

G{z) = I dX^ (16) 

to the upper and lower complex plane, 

G^{z^ X)=Trpjj{X)Ti-^P^{X) (17) 


r[P,m;p] = 

+ J dXp{X) - aln(A) - 7V/ln(A2 -h m^)^ 

J dAdAV(A)p(A')ln(A2 - A'2) 

- - l) j dXplnp (22) 

The P contribution is the self Coulomb subtraction and is 
consistent with the subtraction in the Hilbert transform. 
The saddle point equation 5T/5p = 0 yields the hydro¬ 
static equation ED using the symmetric density p^. 

6. Hydro dynamical instanton. Following the ini¬ 
tial observation in Q, we note that the fixed time zero 
energy solution to m is an instanton with imaginary ve¬ 
locity V = —iA. The conserved current j = pv satisfies 
(r = it) 


so that ED now reads 

aPNz^G - l-PzG^ + zd,G -{a + 27V/)G - 2aPN^z = 0 

(18) 

Defining G{z) = —{2/P)dzlnf{z) and setting z^ = —w 
and f{^/—w) = g{w), we have 


drp{X) = dx [aPNXp{X) - Pttp{X)[ph{X) - ph{-X)] 


a 


+dxp{X) - -p(A) - ^2 


27V/A 


■p{X) 


(23) 


We identify r with the stochastic time, and (1^ describes 
the stochastic relaxation of the chiral eigenvalue density 
of the chiral Dirac spectrum (out of equilibrium) to its 
asymptotic (in equilibrium) hydro-static solution. For 
m = 0, (|23)) can be rewritten as 


2wdyjwg + {—aPNw -I- 1 — a — 27V/) d^g — ^aP^N'^g = 0 

(19) 

We note that ED is Laguerre-like except for the wrong 
sign in the last contribution. The general solution is a 
hyper-geometric function. The spectral density follows 
from the discontinuities of the solution using ED- 


drP^ = aPNdx (Xp^) - Pirdx (p^Ph) 


+dxxp^ 


{a + 2Nf)dx 



After multiplying this equation by A^, we can 
Hilbert transform. The result is 


(24) 


take its 


5. Dyson Coulomb gas. We note that (© can be 
re-written in terms of 1-dimensional Dyson Coulomb gas. 
At large TV the ensemble described by (HD is sufficiently 
dense to allow the change in the measure, 

N 

WdX^zze^'^P^Dp (20) 

i=l 

with iSJA] = f dXp(X)ln(po/p(X)) the Boltzmann en¬ 
tropy [i3- Thus 

Zp[m]^ J ( 21 ) 


drPH = aPNdx {Xpjj) - ^dx 

+dxxPH - (a + 2Nf)dx (25) 

With the usual definition, this gives an equation for the 
time dependent Green’s function 

drG = aPNdz (zG) - ^dz (G2) 

+dzzG -{a + 2Nf)dz Qg^ (26) 

We note that the stationary solution fulfills 

aPNzG-^G'^ + dzG-{a + 2Nf)-G + G = 0 (27) 


with the effective action 
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The constant C is fixed by noting that for large N and 
large z, G/N k. "Ijz so that C = 4aiV^, in agreement 
with (flkll . 

The general time-dependent solution can be analyzed 
by first re-scaling r —>■ t jN. In the large N limit G solves 


X {{dtvf - 


a^G-t;9(^G-az)a^G-a/3G = 0 (28) 

which is a Burgers like nonlinear PDE. We can solve it 
with the method of complex characteristics [l^. For that 
we introduce curves in the space of z and t, parametrized 
by s and labeled by zg, along which the nonlinear PDE 
follows from ordinary ODE 


ds 








(29) 


with the condition that z(t = 0) = zg and t(s = 0) = 0 
(which means t = s). These ODE can be solved for a 
specific initial condition. 


7. Sound waves. To gain some insights to the gen¬ 
eral time-dependent solutions, we analyze first the hydro- 
dynamical equations in the linearized density approxima¬ 
tion. For that, it is convenient to identify the classical 
hydrodynamical action associated to ([5]). Using the stan¬ 
dard canonical procedure we found 


The speed of sound Vs{X) = 2Trj3po{X)/N is local in the 
chiral spectrum. It is Us(0) = 2j5y/a at zero virtuality. 
We note that (l34l) is extensive with N for pq{X)/N nor¬ 
malized to 1. 

8. Chiral relaxation time. An interesting dynam¬ 
ical question regarding the chiral Dirac spectrum is the 
typical relaxation time associated to the formation or dis¬ 
appearance of the spontaneous breaking of chiral symme¬ 
try. This issue is important for QCD with matter under¬ 
going dynamical chiral symmetry restoration or breaking 
close to the chiral transition temperature. Since the set¬ 
up involves off-equilibrium QCD in a moderately strong 
coupling regime, first principle calculations are usually 
challenging. 

Here we answer this question by focusing on the time 
it takes for the eigenvalues near zero-virtuality viewed as 
a hydrodynamical fluid to re-arrange stochastically. For 
that, consider that at r = 0 all the eigenvalues are local¬ 
ized at zero virtuality to mock up an initial and chirally 
restored phase. The relaxation of this phase to a sponta¬ 
neously broken chiral phase is characterized by the time 
it takes the sound wave to fill up the gap or the so-called 
zero-mode-zone (ZMZ) spanned by Wigner semi-circle. 

To show this we use the initial condition G(z = zg, r = 
0) = 2NIzq in (1291) to obtain the solution 


S = y dtdXp{X) {y^ — u[p]) 
which is linearized by 


(30) 


p pz Pq{X)+ 2dx'p and pv pz- dt(p (31) 

Inserting (IHTl) into S yields in the quadratic approxima¬ 
tion 


S 2 = y {{dtV? - pl{X)W'^[ip\) (32) 

with the potential 


W[p] = jdirdx [ph{X) + Ph{-X)] - dx ) ( 33 ) 


G = e“ 


.2N 

zo 


(35) 


The remaining equation is now 


which yields 


dr 


1 


2 

zo 


(36) 


Z = -6“'^" + Uo - 


1 

azg \ azg ^ 

Solving for zg and inserting the answer in (IMI) yield 


(37) 


= ^ (38) 

X ^z\/a± {zxfaY — 4 (l — 


For p{X) = pg(A), after the rescaling Nt —>• t, (1^ sim- The spectral density follows from the discontinuity of 
plifies for large N (fT^ 
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p(A) = ^ (l - ^ (39) 

X (l - 

which is seen to interpolate between a delta-function 
N6{X) at r = 0 and a Wigner semi-circle at asymptotic 
T at a rate given by the speed of sound z;s(0) at zero 
virtuality. The typical time in physical units is 



a _ 7rp(0) _ I (q^q) \ 
Vs{0) ~ 2j3N ~ 2/?n 


(40) 


Recall that in chiral random matrix theory the scale a 
is related to the chiral condensate by the Banks-Casher 
formula 14 ( 9 ^ 9 ) = ~ N^/a. Thus, the last 

equality with n = N/Vi. Near zero virtuality and close 
to the ergodic regime this time is universal. It is non- 
perturbative and fully gauge-invariant. 

In the QCD vacuum we may identify n with the in- 
stanton density [I| (and references therein). Typically, 

I (gig) I Ki 1/(1 fin^) and n k. 1/(1 fm^) so that At k. 
lfm/4 for QCD with Nc = 3 {/) = 2), which is short. 
In matter, both the chiral condensate and the instanton 
density change. For a dyon liquid model in the confin¬ 
ing phase in the range 0.5 < T/T^ < 1, we may identify 
n = ud/Nc with the dyon density [l^. We recall 
that the instanton splits to dyons with only the KK 
dyon carrying the zero-mode. From the analysis pre¬ 
sented in we have for A/. = 2 (/3 = 1) and Nj = 1 


I (qh) I _ / e-^ \ 

2/3nD/fVc ^ T [ aj{T) J 


(41) 


with as{T) = (10/3) ln(T/0.36Tc) the running coupling, 
and C a constant of order 1 in the dyon density which is 
fixed by lattice measurements. 


9. Conclusions. The hydrodynamical description of 
the chiral Dirac spectrum captures some key aspects of 
the stochastic relaxation of the Dirac eigenvalues in the 
diffusive regime. The hydrodynamical set-up supports 
an instanton that describes the stochastic relaxation of 
the Dirac eigenvalues as a fluid. The fluid exhibits sound 
waves that can be used to estimate the time it takes 
for a localized density with chiral symmetry restored, to 
de-localize to a Wigner semi-circle with chiral symmetry 
restored. The relaxation time of a fluid of Dirac eigenval¬ 
ues near zero virtuality captures a typical equilibration 
time for chiral symmetry breaking or restoration in QCD 
that is robust and gauge-independent. 
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